Energy barriers for vortex nucleation in dipolar condensates 
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We consider singly-quantized vortex states in a condensate of ^^Gr atoms in a pancake trap. We 
obtain tlie vortex solutions by numerically solving the Gross- Pitaevskii equation in the rotating frame 
with no further approximations. The behavior of the condensate is studied under three different 
situations concerning the interactions: only s-wave, s-wave plus dipolar and only dipolar. The 
energy barrier for the nucleation of a vortex is calculated as a function of the vortex displacement 
from the rotation axis in the three cases. These results are compared to those obtained for contact 
interaction condensates in the Thomas-Fermi approximation, and to a pseudo-analytical model, 
showing this latter a very good agreement with the numerical calculation. 

PACS numbers: 03.75.Lm, 03.75.Hh, 03.75.Nt 



; I. INTRODUCTION 

Chromium condensates were first experimentally re- 
alized in 2005 1]. In contrast to condensates of al- 
kali atoms, ^^Cr condensates present an additional, non- 
negligible interaction, the dipolar interaction. It is non- 
local, anisotropic and long-range and these features intro- 
duce novel aspects in the physics of Bose-Einstein con- 

■ densates (BECs) 0, Q. The atom-atom interaction is 
then determined by the balance of the s-wave contact 
interaction and the dipolar interaction, allowing thus to 
investigate the effects of their interplay. Although alkali 
atoms also possess some magnetic dipole moment, it is 
thirty-six times weaker than in chromium, which makes 
the latter a suitable system to experimentally study dipo- 

\ lar condensates of atomic gases |l|, 01 • 

The evidence of the superfluid character of BECs 
comes mainly from the appearance of singly-quantized 
vortices above a critical rotation frequency as a means to 
convey angular momentum to the system. Vortices have 
not been experimentally realized in dipolar BECs yet, but 

] they are a most appealing phenomenon and there is in- 
tensive theoretical research on them 0, 0, d, H, [13, [HI ■ 
The critical frequency for vortex nucleation in dipolar 

' condensates has been predicted theoretically in different 
papers [l^, [IH. This critical frequency is based on a 
thermodynamical consideration that takes into account 
the energy of the vortex at the center of the trap p^ . 
However, in the present understanding of vortex forma- 
tion (see Ref. [Tj] for a review), the nucleation process 
takes place at the surface of the condensate and then 
the vortex enters it. Therefore, the theoretical predic- 



tion provides a lower bound to the experimental critical 
frequency 13, . This difference can be understood in 
terms of an energy barrier the system needs to overcome 
to bring the vortex from the surface to the center of the 
condensate [l6l.[T7j. That is, the nucleation of a vortex is 
associated with the existence of an energy barrier in the 
configuration space between the initial vortex-free state 
and the final centered vortex state. The formation energy 
of a vortex can be estimated by calculating the energy of 
a single off-center vortex as a function of the vortex core 
position. 

The aim of this work is to quantify this extra energy 
that should be given to dipolar condensates to nucleate 
a vortex. Since we are interested mainly on the effects of 
the dipolar interaction, we will consider a very small s- 
wave scattering length. In the numerical calculation, we 
solve the three-dimensional (3D) Gross-Pitaevskii equa- 
tion in the rotating frame using the imaginary-time evo- 
lution method. We obtain the vortex states, the vortex 
formation energy and hence the nucleation barrier. 

This work is organized as follows. In SecHIlwe describe 
the theoretical framework and the system under study. 
In Sec. Illll the analytic expressions for the energy barrier 
in the Thomas-Fermi regime are given in 2D and 3D. In 
Sec. lIVI we present the numerical calculation of the nucle- 
ation barrier and compare it both to the Thomas-Fermi 
results for purely s-wave condensates and to a pseudo- 
analytical model which is shown to reproduce very well 
the numerical results. Finally, a summary and conclud- 
ing remarks are offered in Sec. |Vl 



II. THEORETICAL FRAMEWORK 
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Weakly interacting dipolar condensates are well de- 
scribed in the mean-field regime by the Gross-Pitaevskii 
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(GP) equation, which includes a new term taking into 
account the dipole-dipole interaction. In this framework 
it is possible to numerically study yortex states by work- 
ing in the rotating frame [l3|, that is a reference frame 
that rotates at a frequency O around the z axis. The GP 
equation then reads 

- 1 

-^^^ + 14rap(r) + <?|^(r)|2 + l/dip(r) - ^L, i^(r) 
2ra J 



/iV(r) 



(1) 



where m is the atomic mass, ^/'(r) is the condensate wave 
function normalized to the total number of particles, 
is the angular momentum operator along the z axis and /i 
is the chemical potential. The terms in the left-hand side 
of the equation are, respectiyely, the kinetic energy, the 
trapping potential, the s-waye contact interaction, the 
dipolar potential and the rotating potential. When the 
dipolar BEC is at rest, the yortex- free state is obtained 
from Eq. ([T]) by setting = 0. 

The trapping potential is assumed to be harmonic and 
with axial symmetry, being z the symmetry axis 



ap(i") = y (^J 
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where 



+ and uj±^ and lOz 



are the radial and 



axial angular trap frequencies, respectively. 

The contact interaction potential is characterized by 
the coupling constant g — Airll^a/m, with a the s-wave 
scattering length. Using the physics of Feshbach reso- 
nances, the value of a can be tuned to very small values 
by changing the magnitude of an external magnetic field 
[18|. If this tuning is performed far enough from the 
actual resonance, the condensate is not destroyed and 
experiments can be indeed carried out [Toj . 

The mean-field dipolar interaction Vdip(r) is given by 
the integral 



where 



Vdipir) = / dr'i;dip(r-r')|V(r')l' 



fiQfi'^ 1 — 3 cos^ I 
47r |r-r'P 



Udip(r - r') 



(3) 
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is the dipole-dipole potential. In this equation fi is the 
magnetic dipole moment of the atoms, /io is the vacuum 
permeability, r — r' is the distance between the dipoles, 
and is the angle between the vector r — r' and the 
dipole axis, which we also take to be z. Note that dipo- 
lar interactions are anisotropic, being attractive when the 
dipoles are head-to-tail and repulsive when they are side- 
by-side. Therefore, in our symmetry configuration, the 
dipolar interaction will be attractive along z and repul- 
sive in X and y directions. Since we are dealing with 
pancake-like condensates, the overall interaction will be 
mainly repulsive, although the small attractive part is 
what can b ring instability and collapse to the conden- 
sate [23,[2l|,[2l[2l,[ll. 



When the rotation frequency ^ exceeds a critical fre- 
quency flc, the centered-vortex state becomes energeti- 
cally favorable and Eq. (P) converges to such a solution, 
irrespective of the initial conditions. However, an off-axis 
vortex is not a stable state of the system and therefore 
cannot be found as a solution of GP equation as it ap- 
pears in II]). One thus needs to impose some constraint 
during the minimization process to fix the vorticity along 
a line {xy, yv, z). In this work this is achieved as follows. 
First, we generalize the Feynmann-Onsager ansatz for a 
centered-vortex state to the case of an off-axis vortex [il| , 



Xv + i{y - yv) 



^y{x - Xy)'^ + [y-yvY 



(5) 



and use it as the initial wave function. Here 'ipo{r) is the 
ground state vortex-free wave function obtained directly 
from Eq. ([1]) with 17 = 0. This expression fixes a unit of 
vorticity along the line {xy,yy,z) and corresponds to a 
straight vortex line. To maintain such a vorticity during 
the i mag inary-time minimization we impose the condi- 
tions [nj 



Re[V'(a;i,,y, z)] = Vy, z 
Im[?/'(x, ^t,, z)] = Vx, z 



(6) 
(7) 



to the wave function. These relations are the equivalent 
of considering the vortex line as the intersection of two 
fiat nodal surfaces, instead of allowing them to have some 
curvature. With this method, the quantization of the 
circulation is ensured in all cases, but the solutions are 
restricted to the case of straight vortex lines. 

In order to compute the energy of the condensate we 
evaluate the energy density functional in the rotating 
frame, which has the standard GP form but with a new 
term, i^dip, which is the interaction energy due to the 
dipole-dipole potential. 



E[ll}\ — i?kin + ^'trap + -E^int + E, 



dip 



Er = 



(8) 



In this work we will consider a ^^Cr BEC of = 10'°^ 
atoms in a trap with frequencies uj± — 8.4 x 27r s^^ and 
ujz — 98.5 X 27r s~^, with a magnetic dipole moment 
of fi — 6 hb and a scattering length a — S ag. The 
oscillator length in the plane, a± — ^JhJmiD\_^ gives an 
order of magnitude of the size of the condensate, and 
in our case takes the value aj_ = 4.8 /im. By evolving 
the GP equation ([1]) in imaginary-time until convergence 
is reached, we find its solution. This is equivalent to a 
minimization of the energy density functional ([8]). 

In Fig. [T] we show the density profiles of the ground 
state and the singly-quantized vortex state of such a con- 
densate in the three situations: when only s-wave inter- 
actions are considered, when s-wave plus dipolar interac- 
tions are considered, and when only dipolar interactions 
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FIG. 1: Density profiles for vortex (dasfied) and ground 
(solid) states of condensates whose atoms interact via s-wave 
(left panel), s-wave plus dipolar (middle panel) and dipolar 
(right panel) interactions. 



are considered. We see there that the dipolar interactions 
make the condensate more repulsive, hence increasing its 
size in x and y directions. For more detailed discussions 
of condensate deformation due to anisotropic dipolar in- 
teractions see Refs.llJ, 0, Hi and references therein. 



III. THOMAS-FERMI MODELS FOR THE 
ENERGY BARRIER 

In purely s-wave condensates, when the number of par- 
ticles is large and the repulsive contact interaction is 
strong, surface effects at the boundary of the conden- 
sate become negligible. Under these conditions one can 
neglect the kinetic energy term in Eq. ^ and show that 
the density has a parabolic profile. This allows to find 
analytical solutions for the BEC. This approximation is 
known as Thomas-Fermi approximation (TF) and can be 
allied whenever the condition Na/a± ^ 1 is verified 

For dipolar gases there also exists a TF approximation 
[27| . and it is possible to obtain semi-analytical expres- 
sions, which come in terms of a complicated transcen- 
dental relation. In this case the condition Na/a± ^ 1 
does not define the TF regime, since it does not take the 
dipolar interaction into account. However, a qualitative 
characterization of the TF limit in dipolar condensates 
relies on the parabolic density profile of the ground state 

M- 

For the sake of clarity, we will specify when we are con- 
sidering TF results only applicable to s-wave condensates 
(SWTF), or when we are using a more general approach 
which only relies on the parabolic profile of the conden- 
sate and whose results can be applied both to s-wave and 
dipolar cases (TF). 

The energy barrier AE{d, ft) that has to be overcome 



to nucleate a vortex is given by the vortex energy in the 
rotating frame as a function of the vortex distance d = 
y/xl+y^ from the symmetry axis. It can be calculated 
from the energy in the laboratory frame E{d) using 



AE{d,n) = E{d) - nL,{d) 



(9) 



where is the expectation value of the angular momen- 
tum operator L^. To simplify the notation, we consider 
AE{d, n) and E{d) to be the energies referred to the 
ground state of the system. Note that the energy bar- 
rier depends on the actual value of but the energy in 
the laboratory frame does not. For = ^lc and d = 
we have A£'(0, ilc) = since under this conditions the 
energy of a centered vortex is exactly the rotation en- 
ergy, which comes directly from the definition of critical 
frequency [l^ . For a vortex near the boundary of the 
condensate ^ [1^, whence AE{d, fl) 0. Between 
these two limits, the energy barrier reaches a maximum 
value AE^max at the position dmax, corresponding to the 
extra energy the system needs in order to carry the vortex 
from the surface to the center of the condensate. 

In the TF limit it is possible to obtain analytic expres- 
sions for the energy barrier for vortex nucleation. Al- 
though we are dealing with 3D systems, since the trap 
anisotropy parameter is large, A = — 11, it is 

worth considering the 2D limit. We will derive an ana- 
lytical expression for 2D systems and briefly discuss the 
results for 3D condensates [iJ, [iQ, [13 . 



A. Two-dimensional case 

In the TF limit the vortex core can be considered to 
be much smaller than the radius of the condensate, i.e. 
^ ^ R±, where ^ is the healing length. Then, a good 
approximation is to consider the energy of the vortex 
to be proportional to the energy of the centered-vortex 
configuration [s^l 



Pjd) 
P(0) 



m 



(10) 



where p{d) and p(0) are, respectively, the ground state 
densities at positions d and d = 0. Then the energy 
barrier takes the form 



AE{d,n) 



piO) 



E{0) - nL,{d) 



(11) 



To obtain an expression of the barrier, one substitutes 
p(d), p{0) and Lz{d) in Eq. (fTT|) by their TF expressions 
[29|, which assume a parabolic profile of the density of 
the vortex-free state. This procedure yields 



AE-'^id.Q) 



2D 



Rtf 



\Rtf 



(12) 
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Maximizing Eq. ([T2)) . one can find the position of the 
barrier height 



- Rtf\I 1 - 



2n 



(13) 



and its value 

AEl^^^in) ^ AE'^idn..., n) = , (14) 

where the TF radius in the radial direction Rrf is related 
to its root-mean-square (rms) analogue R± by Rtf = 
V5R±. 

In the SWTF limit, one can find an analytical expres- 
sion for the critical frequency [3l| 



2D 



mR'^p 



In 



O.SSSRtf 



(15) 



which takes the value fl^^ — 0.34 wj^ for our parameters. 
To obtain Eq. (fTS]) an explicit expression for ^ is used 
which can only be applied to s-wave condensates. 



B. Three-dimensional case 

In 3D, the procedure to obtain an expression for the 
energy barrier differs from the 2D case. It has been re- 
ported elsewhere [l^, [l^ [l3l , and here we only summarize 
the main results. The nucleation energy barrier has the 
expression 



d 



3D 



1 



Rtf 



3/2 



Rtf 



5/2 



with the maximum at the position 



c^m^x - Rtf\I 1 - 



3r!P 



5 n 



(16) 



(17) 



and a barrier height 



(18) 

As in the 2D case, the relation Rtf — V^R± also holds. 

Note that expression ([16]) has been found only using 
that the density profile is parabolic, and makes no specific 
assumption of purely s-wave interacting BECs. There- 
fore, it also applies when dipolar interactions are consid- 
ered, provided the density profile is parabolic, as in 2D, 

Eq. ini). 

Again, in the SWTF limit one can find an analytic 
expression for the critical frequency [3l| 



B n^jmiRrF 



2 mR^p 



(19) 



which takes the value — 0.38 lo± m the present work. 
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FIG. 2: (Color online) Energy barrier for the nucleation of 
a vortex in the rotating frame at f2 = fie as a function of 
the vortex displacement from the center. The dashed line 
corresponds to the pure contact interaction BEC (with a — 
Sag), the solid line corresponds to a pure dipolar BEC, and 
the dash-dotted hne to a condensate with contact plus dipolar 
interactions. 



IV. NUMERICAL CALCULATION OF THE 
ENERGY BARRIER 

In this section we use the numerical procedure outlined 
in Sec. |TT] to solve the GP equation for different vortex 
distances d from the symmetry axis. An estimate of the 
formation energy of the vortex is then obtained from the 
energy in the rotating frame given by the energy func- 
tional Eq. ([8]) relative to the ground state energy. 

We plot in Fig. [2] the vortex formation energy as a 
function of the vortex displacement from the center, for 
the same cases as shown in Fig. [T]and for 17 = fie- The 
critical rotation frequency is calculated from thermody- 
namical arguments jl2l . [l3| , giving 



1 



NL. 



-{E{Q)-Eo) , 



(20) 



where E{0) and Eq are the energies of the centered vortex 
state and the ground state in the laboratory frame, re- 
spectively. We see in the figure that the centered- vortex 
state and the vortex-free state have the same energy in 
the rotating frame but they are indeed separated by an 
energy barrier. The critical frequencies are different in 
each case, namely: flc = 0.45?ia-'_L for a purely s-wave 
condensate, and ilc = 0.25huj± for a condensate with 
both interactions and for a purely dipolar condensate. 

The dashed line in Fig. [5] corresponds to the purely 
contact interaction BEC, the dash-dotted line to a con- 
densate with contact plus dipolar interactions, and the 
solid line corresponds to a purely dipolar BEC. The effect 
of the strength of the overall interaction is to decrease the 
barrier height: it is higher for an s-wave condensate be- 
cause it is the weakest interacting case (see Fig. [T|) and 
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s-wave 

both 

dipolar 

3D SWTF 
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TABLE I: Comparison of the energy barrier maximum and 
its position obtained in the numerical simulation (sim.) and 
the 2D and 3D pseudo- analytical results. 
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FIG. 3: (Color online) Energy barrier for the nucleation of a 
vortex in the rotating frame at f2 = fJc for the same cases as 
in Fig. (2] together with the predictions of the SWTF model 
in 2D and 3D. 



smaller for a condensate with dipolar plus s-wave inter- 
actions, since here the effective repulsion is the strongest. 
Besides, when dipolar interactions are considered the ra- 
dial extent of the condensate is larger and thus the barrier 
height is displaced towards the surface, as compared to 
the purely s-wave case. 

To get rid of the effect of the different sizes due to the 
different strengths of the repulsive interactions, we plot 
the energy barriers in Fig. [3] as a function of the vortex 
displacement normalized to the transversal rms-radius of 
the corresponding ground state, d/R±. It is interesting to 
note that even though R± is different for each curve, the 
three critical barriers have the same qualitative behavior 
when expressed in units of R± . 

In Fig. [3] we compare the numerical barriers to those 
obtained from the SWTF results, which correspond to 
Eqs. (fT2|) and (fT6|) using the critical frequencies Eqs. p5|) 
and Uni), respectively. The disagreement between the 
SWTF prediction and the numerical calculations can be 
understood recalling that the parameters used in the 
present work lead to Na/a± — 5.5, which is not much 
larger than 1 and therefore we are not in the full range 
of validity of the SWTF approximation. 

Nonetheless, to obtain an order of magnitude for the 
energy barriers, it is possible to combine the TF results 
(fT2|) and (|16p with the numerical values of fie and R± 
obtained from the GP equation. This pseudo-analytical 
method avoids thus the calculation of the complete en- 
ergy barrier. Table |T] shows the main parameters of the 
barriers using these approximations at the critical fre- 
quency. 

We see from Table U that the 3D approximation is in 
good agreement with the characteristics of the barriers 
found in the numerical simulation, whereas the 2D ap- 
proximation shows some deviation. This means that in 
spite of the large anisotropy of the trapping potential, 
the 3D nature of the condensate is important to charac- 
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0.25 
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2D 
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1.22 








3D 
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1.09 



terize the vortex energetics. In addition, the energy of 
the maximum is better reproduced for condensates with 
dipolar or both interactions. This happens because the 
corresponding density profiles (see Fig. [1]) are closer to 
the TF inverted parabola profile than those of the purely 
s-wave condensate. 



V. CONCLUSIONS 

In this work we have studied the energy barriers for 
vortex nucleation in the case of dipolar condensates. We 
have considered the regime of small scattering lengths, 
which in our case means that dipolar effects are en- 
hanced, and have obtained the vortex formation energies 
as a function of the vortex displacement from the sym- 
metry axis. We have compared the results obtained when 
s-wave plus dipolar interactions are considered with those 
obtained in the cases of condensates interacting only via 
s-wave interactions and condensates with only dipolar 
interactions. 

Completely analytic expressions can be found in the 
s-wave Thomas-Fermi regime for the energy barriers in 
2D and 3D. From these expressions the position of the 
barrier maximum and its value can be obtained. Since 
they mainly hold for Na/a± 3> 1, their predictions are 
not in much agreement with our numerical results. 

However, a pseudo-analytical approach can be used 
which is based on the equations from TF theory, but 
taking the critical frequency and radius of the condensate 
from the numerical results. This procedure only assumes 
parabolic density profiles for the ground state densities, 
thus being applicable to both s-wave and dipolar conden- 
sates, provided they satisfy the above condition. 

For condensates with dipolar interaction, this pseudo- 
analytical method has been shown to provide close quan- 
titative results for the nucleation energy barrier when 
the 3D Thomas-Fermi expression is used. Since the den- 
sity profiles for the dipolar and s-wave plus dipolar con- 
densates are approximately parabolic, the Thomas-Fermi 
expression for the energy barrier holds and its predic- 
tions are in accordance with the numerical results. This 
method can be a useful tool to estimate the energy bar- 
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rier for vortex nucleation, avoiding thus the calculation 
of the complete energy barrier. 
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